Damping of Quantum Vibrations Revealed in Deep Sub-barrier Fusion 
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We demonstrate that when two colliding nuclei approach each other, their quantum vibrations are 
damped near the touching point. We show that this damping is responsible for the fusion hindrance 
phenomena measured in the deep sub-barrier fusion reactions. To show those, we for the first time apply 
the random-phase-approximation (RPA) method to the two-body '^O + '^O and ^''Ca + "^''Ca systems. 
We calculate the octupole transition strengths for the two nuclei adiabatically approaching each other. 
The calculated transition strength drastically decreases near the touching point, strongly suggesting the 
vanishing of the quantum couplings between the relative motion and the vibrational intrinsic degrees of 
freedom of each nucleus. Based on this picture, we also calculate the fusion cross section for the '^''Ca 
+ ^°Ca system using the coupled-channel method with the damping factor simulating the vanishing of 
the couplings. The calculated results reproduce well the experimental data, indicating that the smooth 
transition from the sudden to adiabatic processes indeed occurs in the deep sub-barrier fusion reactions. 
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Heavy-ion fusion reactions at low incident energies serve 
as an important probe for investigating the fundamental 
properties of the potential tunneling of many-body quan- 
tum systems. When two nuclei fuse, a potential barrier, 
called the Coulomb barrier, is formed because of the strong 
cancellations between the Coulomb repulsion and the at- 
tractive nuclear force. The potential tunneling at inci- 
dent energies below this Coulomb barrier is called the sub- 
barrier fusion. One important aspect of the sub-barrier fu- 
sion reactions is couplings between the relative motion of 
the colliding nuclei and nuclear intrinsic degrees of free- 
dom, such as collective vibrations of the target and/or pro- 
jectile lUl]. Those couplings result in the large enhancement 
of the fusion cross sections at the sub-barrier incident en- 
ergies, as compared to the estimation of a simple potential 
model. The coupled-channel (CC) model taking into ac- 
count the couplings has been successful in explaining this 
enhancement L2>^. 

Recently, it has been possible to measure the fusion cross 
sections down to extremely deep sub-barrier incident en- 
ergies |0-|3l- The unexpected steep fall-off of the fusion 
cross sections, compared to the standard CC calculations, 
emerges at the deep sub-barrier incident energies in a wide 
range of mass systems. These steep fall-off phenomena are 
often called the fusion hindrance. The emergence of the 
fusion hindrance shows the threshold behavior, which is 
strongly correlated with the energy at the touching point of 
the two colliding nuclei In this respect, it has been 

shown that a key point to understand this fusion hindrance 
is the potential tunneling in the density overlap region of 
the two colliding nuclei (see Fig. 1 in Ref. [8]). 

To describe the fusion hindrance phenomena, many the- 
oretical models to extend the standard CC model have been 
proposed. Based on the sudden picture, Mi§icu and Es- 
bensen have proposed that a strong repulsive core exists in 
the inner part of the Coulomb barrier due to nuclear incom- 



pressibility ifToll . This model can reproduce well the fusion 
hindrance from the light- to heavy-mass systems [lOl - 
fT2|l. Dasgupta et al. proposed the concept of the quan- 
tum decoherence of the channel coupling s ol , but there are 
only simple calculations with this model |13j. Based on the 
adiabatic picture, which is the opposite limit to the sudden 
approach, Ichikawa et al. introduced the damping factor in 
the standard CC calculations in order to smoothly joint be- 
tween the sudden and adiabatic processes [14.1 . This model 
can reproduce the fusion hindrance better than the sudden 
model. However, the physical origin of the damping factor 
was still unclear. 

In this Letter, we show the physical origin of the damp- 
ing factor proposed in Ref. [|l4ll . In the standard CC model, 
it has been assumed that the properties of the vibrational 
modes do not change, even when two colliding nuclei touch 
with each other. However, in fact, the single-particle wave 
functions drastically change in the two nuclei approach- 
ing each other. This results in the damping of the vibra- 
tional excitations, that is, the vanishing of the couplings 
between the relative motion and the vibrational excitations 
of each nucleus. To show this, we for the first time ap- 
ply the random-phase approximation (RPA) method to the 
two-body '^O -i- '^O and '^''Ca -i- "^°Ca systems and calculate 
the octupole transition strength, B(E3), as a function of the 
distance between the two nuclei. We below show that the 
obtained B(E3) values for the individual nuclei are indeed 
damped near the touching point. 

To illustrate our main idea, we first discuss a disappear- 



ance of the octupole vibration during the O -i- O fu- 
sion process. We calculate the mean-field potential with 
the folding procedure using the single Yukawa function to 
conserve its inner volume [15]. In the two-body system be- 
fore the touching point, we assume the two sharp-surface 
spherical nuclei. After the touching point, we describe the 
nuclear shapes with the lemniscatoids parametrization, as 
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FIG. 1. (color online) Nilsson diagram for the neutron single- 
particle states versus the distance between "'O + '^0. The solid 
and dashed lines denote the positive and negative party states, 
respectively. The gray area denotes the overlap region of '^O + 
'^O. The solid and dashed arrows denote the main p-h excitations 
generating the octupole vibration. The single-particle density dis- 
tributions for these p and h states are given in the insets from (a) 
to (d) and (a') to (d'). 



shown in Ref. fl^. Using this, we can describe the smooth 
transition from the two- to one -body mean-field potentials. 
The depths of the neutron and proton potentials are taken 
from Ref. [17]. We use the radius for the proton and neu- 
tron potentials, Rq, with Rq = 1.26/l'^^ fm, where A is the 
total nucleon number. In the calculations, the origin is lo- 
cated at the center-of-mass position of the two nuclei. The 
above procedure works well, because Umar and Oberacker 
found that the double -folding potential with the frozen den- 
sity agrees almost perfectly with the density-constrained 
time-dependent Hartree-Fock (TDHF) approach for dis- 
tances ^ > 6 fm Hi]. 

Using the obtained mean-field potentials, we solve the 
axially-symmetric Schrodinger equation with the spin- 
orbit force. Then, the parity, jt, and the z component of 
the total angular momentum, Q, are the good quantum 
numbers. The details of the model and the parameters are 
similar to Refs. lIlSl IITII . We calculate the single-particle 
wave functions of both the projectile and target using the 
one-center single Slater determinant. We expand the total 
single-particle wave function by many deformed harmonic- 
oscillator bases in the cylindrical coordinate representation. 
The deformation parameter of the basis functions is deter- 
mined so as to cover the target and projectile. The basis 
functions are taken with its energy lower than 25 Jico. 

Figure 1 shows the Nilsson diagram for the obtained neu- 




FIG. 2. (color online) Transition densities and currents for the 
first excited 3" state with Q." = 0^ atR = (a) 14.0 fm, (b) 8.0 fm, 
and (c) 6.4 fm. The contour lines denote the transition density. 
The arrows denote the current density. These two values are nor- 
malized in each plot. The (red) thick solid line denotes the half 
depth of the mean-field potential. 



tron single-particle energies versus the distance between 
16q _|_ 16q rpjjg solid and dashed lines denote the positive 
and negative party states, respectively. The gray area de- 
notes the overlap region of the two nuclei. The distance 
R = 6.4 fm corresponds to the touching point. Some den- 
sities for the obtained single particles at i? = 14 and 6.4 
fm are given in the insets from (a) to (d) and (a') to (d') in 
Fig. 1, respectively. AtR = 14 fm, we see that the positive- 
(negative-) parity indicates the symmetric (asymmetric) 
combinations of the single-particle states referring to the 
right- and left-sided '^O. Thus, the positive- and negative- 
parity single-particle states are degenerate for large R. With 
decreasing R, these single-particle states smoothly change 
to those for the composite ^^S system. 

We can now easily extend the RPA method [19] to the 
two-body system, because the wave functions of both the 
one- and two-body systems are described with the one- 
center Slater determinant. We can directly superpose all 
combinations of the particle (p) and hole (h) states for the 
obtained single particles in a unified manner for both the 
one- and two-body systems. We solve the RPA equation at 
each center-of-mass distance between '^O -i- ^^O. At large 
R values, the RPA solutions with fl" = 0^ and 0~ represent 
the symmetric and asymmetric combinations of the states 
where the RPA modes are excited in either the right- or 
left-sided '^O. When R decreases below the touching point, 
they smoothly change to excitation modes in the composite 
^^S system. In the calculations, we only take into account 
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FIG. 3. Transition strengths, B(E3), for the first excited state 
versus the distance between (a) ^^O + '^0 and (b) "^''Ca + '^''Ca. 
The solid line denotes the calculated results. The gray area de- 
notes the overlap region of the two colliding nuclei. In the bottom 
panel (b), the dotted line denotes the damping factor of Eq. (2) in 
Ref. [141] estimated from the experimental data of the fusion cross 
section. 



the p-h states with the excitation energies below 30 MeV. 
We use the residual interaction as the density-dependent 
contact one taken from Ref. [20J. The strength of the resid- 
ual interaction is determined at each R such that the lowest 
D." = 0~ solution of the RPA appears at zero energy. We 
have developed a new scheme based on the Tomonaga the- 
ory of collective motion |2ll 221, that enables us to separate 
the center of mass, the relative motion, and the intrinsic de- 
grees of freedom. It is a generalization of the known pro- 
cedure in the RPA [23]. It is also interesting to compare 
our results with the TDHF dynamics i24ll . 

For the calculated results, the obtained first-excited state 
of the right-sided '^O is the octupole (3~) one with a large 
B(E3) value. At 7? = 15 fm, the excitation energy and the 
BiE3, 3- Op value are 5.29 MeV and 102.07 e^ fm^ 



respectively. We have checked that those values are consis- 
tent with the calculated results of the one- body i^O. Pi gure 
2(a) shows the calculated transition densities and currents 
U% for the first excited state with Q'' = 0"^ at ^ = 14 fm. 
In Fig. 2(a), we can clearly see the octupole vibrations in 
both ^^O's. At i? = 8.0 fm, the transition density of the 
neck part between two '^O's develops [see Fig. 2(b)]. At 
R = 6.4 fm, the octupole vibrations of each '^O become 
weak [see Fig. 2(c)]. The degenerating excitation energies 
of the first-excited states with Q." = 0"^ and 0~ split below 
/? = 8 fm. They become 5.82 and 4.83 MeV at^ = 6.4 fm. 

To more clearly see the damping of the octupole vibra- 
tions, we calculate the B(E3) value for the right-sided '^O. 
We can easily calculate it by taking a symmetric linear 
combination of the octupole transition amplitudes for the 
positive- and negative-parity RPA solutions. Figure 3 (a) 
shows the calculated B(E3, 3^ — » Op values versus the dis- 
tance between '^O + ^^O. In Fig. 3 (a), we can see that the 



B{E3) value for the right-sided '^O (the solid line) falls off 
at around R = S fm with decreasing R, indicating that the 
octupole vibrations are strongly suppressed near the touch- 
ing point. 

The damping of the vibrations originates from the 
change of the single-particle wave functions. At R = 14 
fm, the major p-h excitations generating the octupole vi- 
bration are those from the pi/2 to the d5/2 single-particle 
states in '^O. In these states, the degenerating positive- and 
negative-parity doublet states contribute equally to gener- 
ate the octupole vibration. Those can be seen in the den- 
sity distributions of the fl" = I /2'^ and 1 /2~ states given 
in the insets from (a) to (d) in Fig. 1. The correspond- 
ing p-h excitations are denoted by the solid arrows from 
(a) to (b) and (c) to (d) in Fig. 1. When the two nuclei 
approach each other, the features of these single-particle 
wave functions drastically change. AtR = 6.4 fm, the neck 
formation takes place in the positive-parity states, while it 
is forbidden for the negative-parity states (having nodes at 
the touching point). Thus, the density distributions of those 
parity partners become quite different from each other, [see 
the insets from (a') to (d')]. In the RPA calculation, the 
contributions from the negative-parity states to the octupole 
vibration become small at the touching point [see the solid 
and dotted arrow from (a') to (b') and (c') to (d'), respec- 
tively], resulting in the decreases of the B(E3) value. 

The mechanism for the damping of the quantum vibra- 
tion is a general one valid also for heavier mass systems. 
We have also performed the RPA calculations for the '^''Ca 
+ '^''Ca and ^^Ni + ^^Ni systems. In the calculations, we 
use Rq = l.llA^I^. We obtained the similar damping of the 
B(E3) values for the both systems. In Fig. 3(b), the solid 
line denotes the calculated result for the '^''Ca -i- ^^Ca sys- 
tem. The calculated excitation energy and B(E3, 37 —> Ot) 



value for the '^"Ca +'*"Ca system are 3.24 MeV and 1253.17 
fm^ atR = 15.0 fm, respectively. 

As shown above, the octupole vibrations are damped 
near the touching point, resulting in the vanishing of the 
couplings between the relative motion and the vibrational 
intrinsic degrees of freedoms. This vanishing would lead 
to the smooth transition from the sudden to adiabatic pro- 
cess, as shown in Ref. lll4ll . It is clear that such effect has 
not been taken into account in the standard CC model. One 
candidate to include it is the introducing of the damping 
factor proposed in Ref. [[l4ll . To clearly see the effect of the 
vanishing of the couplings, we calculate the fusion cross 
section for the '^''Ca + ^"Ca system using the computer code 



CCFULL 112511 coupled with the damping factor based on the 
model of Ref. M . 

In the calculations, we include the couplings to only the 
low-lying 3~ and 2"^ states and to single phonon and all mu- 
tual excitations of these states. We take the energies and the 
deformations of each state taken from Ref. [7] to reproduce 
well the experimental data. We use the same deformation 
parameters for the Coulomb and nuclear couplings. For the 
parameters of the Yukawa-plus-exponential (YPE) model, 



4 



Exp. I — • 

No coupling 

YPE 

+ damping 

^/o ■ ! I I I I I I I I I I I I I I 




65 



deed damped near the touching point. We have also shown 
that the damping factor proposed in Ref. lll4j| strongly cor- 
relates with the calculated B(E3) values and the calcula- 
tions of the fusion cross section coupled with the damp- 
ing factor reproduce well the experimental data. This indi- 
cates that the fusion hindrance originates from the damping 
of the quantum couplings and strongly suggests that the 
smooth transition from the sudden to adiabatic processes 
occurs near the touching point. 
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FIG. 4. (color online) Fusion cross sections (upper panel) and its 
logarithmic derivative (lower panel) for the '*''Ca + ^"Ca system 
versus the incident energies. The solid circles denote the experi- 
mental data taken from Ref. The solid and dashed lines de- 
note the calculated results of the coupled-channel method using 
the YPE potential with and without the damping factor, respec- 
tively. The dotted line denotes the calculated result without the 
couplings. 

we use ro = 1.191 fm and a = 0.68 fm. 

It is remarkable that the damping factor strongly corre- 
lates with the calculated B(E3) value for the '^''Ca + '*°Ca 
system. To show this, we take /Iq, = in the damping factor 
of Eq. (2) in Ref. lfl4ll for simplicity. We tune the parame- 
ters Rd and in the damping factor so as to reproduce the 
experimental data of the fusion cross section. We obtain 
Rd = 9.6 fm and = 0.9 fm as the best fit to the data. 
In Fig 3(b), the dotted line denotes the obtained damping 
factor normalized atR= 15 fm. 

Figure 4 shows the calculated fusion cross section (upper 
panel) and its logarithmic derivative dln(Ec.m.o'fus)/dEc.m. 
(lower panel) for the ^"Ca -i- ^''Ca system. The solid and 
dashed lines denote the calculation with and without the 
damping factor, respectively. The dotted line denotes the 
calculation without the couplings. In Fig. 4, we can see 
that the calculated results with the damping factor repro- 
duce well the experimental data, which is better than the 
sudden model JJO- In our model, the calculated astro- 
physical S-factor has a peak structure. We also performed 
the CC calculation for the "^^Ca -i- ^^Ca system and the cal- 
culated result reproduces well the experimental data. The 
CC calculations with the damping factor also already re- 
produced well the experimental data for the ^"^Ni -i- ^^Ni, 

58Ni + 58j^j^ ^jjj 16q _^ 208p|, gygtemS [14]. 

In summary, we have demonstrated the damping of 
the quantum vibrations when two nuclei adiabatically ap- 
proach each other To show this, we for the first time ap- 
plied the RPA method to the two-body '^O + '^O and "^^Ca 
-I- '^''Ca systems and calculated the B(E3) values of each nu- 
cleus. We have shown that the calculated B(E3) value is in- 
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